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Abstract 

The feasibility of shell-model calculations is radically extended by the Quan- 
tum Monte Carlo Diagonalization method with various essential improve- 
ments. The major improvements are made in the sampling for the generation 
of shell- model basis vectors, and in the restoration of symmetries such as an- 
gular momentum and isospin. Consequently the level structure of low-lying 
states can be studied with realistic interactions. After testing this method 
on 24 Mg, we present first results for energy levels and E2 properties of 64 Ge, 

indicating its large and 7-soft deformation. 
PACS numbers: 21.60.Ka, 21.60.Cs, 24.10.Cn 
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The nuclear shell model has been successful in the description of various aspects of 
nuclear structure, partly because it is based on a minimum number of natural assumptions. 
Although the direct diagonalization of the hamiltonian matrix in the full valence-nucleon 
Hilbert space is desired, the dimension of such a space is too large in many cases, preventing 
us from performing the full calculations. The direct diagonalization has been carried out up 
to 48 Cr [jlj. Recently, in order to relax this restriction drastically, stochastic approaches, for 
instance, the shell model Monte Carlo (SMMC) method 0, have been investigated. In fact, 
ground-state || and thermal properties @] have been described well by the SMMC method. 

We have presented the Quantum Monte Carlo Diagonalization (QMCD) method || by 
utilizing the auxiliary field Monte Carlo technique as in the SMMC method, but in a quite 
different way. In general, low-lying states of nuclei are described to a good extent in terms 
of static and/or dynamic mean fields and their fluctuations. The basic idea of the QMCD 
method is to diagonalize the shell model hamiltonian, by using this property, in a subspace 
spanned by a small number of selected basis states obtained by stochastically generated 
one-body fields. Thus, the ground state and several excited states can be obtained. The 
QMCD method has been applied to the Interacting Boson Model In this Letter the 

QMCD method is revised considerably in various aspects so as to be capable of performing 
large-scale shell model calculations with realistic nuclear forces. As examples, 24 Mg and 
64 Ge are taken. In particular, 64 Ge is an N = Z proton-rich unstable nucleus manifesting a 
7- soft structure, with a wide range of theoretical interpretations ( see ref. [|7| ). Thus, the 
shell model calculation can play a crucial role for clarifying the level structure, but so far 
such attempts have been impossible due to the large dimension (~ 1 x 10 9 ). 

We first sketch the QMCD process very briefly, referring to relevant equations of ref. 
0. More details on certain basic points can be found in ref. ||. The shell model hamilto- 
nian consisting of single particle energies and a two-body interaction can be written in the 
quadratic form of Nf one-body operators O a : 

H=J2(E a O a + -V a 2 a ). (1) 
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We consider the imaginary time evolution operator with N t slices: = E[n=i e 5 

where A/3 = (3/N t . By applying the Hubbard- Stratonovich transformation at each time 
slice HH, e _/3H can be expressed as the integral of an operator, n^=Li e~ A/3/l ^ f7 ™\ over Nf x 
N t auxiliary fields a an (see Eq.(4) of ref. j5|) with the Gaussian weight factor G(a) = 
exp (— J2 a ,n A(3/2 | V a I <7an)- The one-body hamiltonian h(a n ) is defined by 



where s a = ±1 (= ±i) if V a < (> 0). In the QMCD method, by generating a new set of 
values for o = {a an } stochastically according to G(a), a new many-body state is created as, 

Nt 

|$(cx)> oc J] e - A ^) |^(o) ); (3) 

n=l 

where |\1/^ ^) is an appropriate initial state. The Hamiltonian is diagonalized in the Hilbert 
subspace spanned by this state and the basis states previously obtained. If this new state 
improves the result of the diagonalization sufficiently well, this state is added to the basis 
states. The number of such basis sates is referred to as the QMCD basis dimension, and is 
increased until reasonable convergence is achieved. 

It is convenient to adopt basis states in the form of Slater determinants: IIa=i a a I — )> 
where N denotes the number of valence nucleons, | — ) is an inert spherical core, and a) a 
represents the nucleon creation operator in a canonical single-particle state a, which is a 
linear combination of the spherical bases. Note that, if | \]>( )) is a Slater determinant, | $(cr)) 
in Eq.(^) remains in the form of a Slater determinant. 

While the QMCD method outlined so far is applicable to fermion systems, its capability 
is limited to simple cases, for instance, a single-j model. Difficulties arise, for example, due 
to finite single particle energies. Thus, a substantial further improvement of the method 
is required for realistic shell model calculations. Such improvements are: (i) the sampling 
scheme is modified, and (ii) additional processes are included to restore symmetries. 

We start with the sampling. In the original version of the QMCD method, a rather naive 
sampling is performed (see Eq.(4) of ref. f|). This sampling creates many unnecessary basis 



(2) 



a 
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vectors in general, and indeed the actual sampling has to be modified for large-scale realistic 
shell model calculations so that important basis vectors are generated still stochastically but 
more efficiently by considering the many-body dynamics. 

The modification regarding the sampling consists of two parts. In the first part, the basis 
state generation is refined so as to make use of the local Hartree-Fock (HF) energy minima. 
In the QMCD calculation, one has to generate good basis states (i.e., Slater determinants 
in deformed bases) which have (i) low values of diagonal matrix elements and/or (ii) large 
off-diagonal matrix elements of the hamiltonian. The point (i) can be fulfilled by using, as 
| \|>( )) in Eq.(^|), a deformed HF solution within the present shell model space. The QMCD 
process is comprised practically of several segments starting with different initial states, 
which are HF states at different local minima. States around a minimum satisfy point (ii) 
in most cases. We then rearrange the one-body evolution process so that the basis states 
are sampled most frequently near the HF local minima, accelerating the generation of state 
vectors having larger overlap with eigenstates of interest. 

The hamiltonian is rewritten, by introducing the constants c a , as 

H = J2(E a O a + \v a {O a - c a ) 2 + V a c a O a ), (4) 

a * 

where a constant term is omitted. After the HS transformation, the one-body hamiltonian 
becomes 

K^n) = J2(( E <* + v » c a)O a + s a V a a an O a ), (5) 

a 

where the c-number — J2 a SaVaVanCa is omitted since it does not change the wave functions 
apart from the normalization. In this expression, the modified one-body term J2 a {E a + 
V a c a )O a includes effects coming from the two-body interaction. The Cq- S £1X6 taken in such 
a way that this term becomes the HF single-particle hamiltonian, Huf- With this /ihf, the 
QMCD basis state takes the form, 

Nt 

|$(cr)) oc Y[ e" A/3(%F+ ^ s <±v a a an o a ) |^,(o)^ ( 6 ) 

n=l 
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Thus, we simply replace the single particle energy J2 a ^aO a by /i H f- If \^^) is the HF 
state being considered, the sampling around a = generates various states around this HF 
state, including Tamm-Dancoff-type states to first order in a, and so on. This treatment is 
possible for all HF local minima. 

In cases of non-spherical nuclei, many Hartree-Fock local minima appear in the search 
for the initial state. By stochastically taking those local minima as the initial states, it is 
possible to take into account a wider variety of configurations. 

The second major improvement on the sampling is the ordering of one-body fields ac- 
cording to their importance. The QMCD method is a method for generating favorable basis 
states for diagonalization, and there is no need to carry out the stochastic integration over 
all auxiliary fields. In constructing the basis states, we start with the most relevant part 
of the hamiltonian, which yields fewer auxiliary fields than the whole hamiltonian. The 
calculation can then be performed more efficiently After certain basis states are obtained, 
we take an enlarged portion of the hamiltonian, so that other terms of the hamiltonian can 
be properly included in constructing the basis states. Eventually the completeness of the 
QMCD basis is guaranteed for the ground state by taking all fields. 

In most cases, the auxiliary fields with large values of | V a | in Eq.(|I]) turn out to have 
quadrupole, hexadecapole or monopole nature. Therefore it is reasonable to arrange all fields 
in descending order of | V a |, and take them starting from the largest one. In addition, for 
a fixed initial state | \&(°)), the total strength of each O a changes due to the Pauli principle 
and to collective effects. Therefore we consider an excitation sum-rule; 

S a = \ OlO a |tf(o))_ |(*(o) \O a | 2 , (7) 

and use it as a practical measure of importance of the O a 's. The selection of O a 's according 
to | V a | and S a plays an essential role in the actual calculations. 

Incorporating all the above improvements, the sampling is made much more efficient. 
Note that this way of sampling clearly differs from that of the SMMC 

We now come to the restoration of symmetries. We implement explicitly kinematic sym- 
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metries such as angular momentum and isospin into the QMCD method, since the restoration 
of such symmetries proceeds only very slowly for wave functions generated stochastically. In 
the previous paper || we have presented the M-projection method to restore the magnetic 
quantum number. 

Since a nucleus has rotational symmetry, the restoration of the total angular momentum, 
denoted as J, is quite crucial. In the QMCD method, we diagonalize the hamiltonian in 
the laboratory frame by using QMCD bases. If the QMCD bases contain all components 
(i.e., Slater determinants) required for the coupling to a good angular momentum, the 
diagonalization restores the rotational symmetry. We accelerate this restoration process, by 
considering rotated states exp(—i9 y J y ) exp(—i9 z J z ) | $(cr)) as candidates of new basis states. 
We have found that the restoration of the angular momentum is remarkably improved by 
taking only several values of the angle #'s. We refer to this method as J-drive. In addition 
to this, the M-projection is carried out for all bases thus created. 

We next discuss isospin. The isospin projection is possible in the same way as the J- 
projection. In this Letter, however, since we consider only N=Z nuclei, we keep good isospin 
in an alternative way. In the decomposition process, Eq. QT]) , all one-body operators can be 
chosen so as to carry a definite isospin T = or 1 for the isoscalar hamiltonian. Since the 
isoscalar fields are dominant over the isovector ones, particularly, for T=0 states, we start 
the QMCD basis generation process with the isoscalar fields. Thus, since the initial HF 
state has T = 0, the isospin is conserved at least until the isovector fields are activated. It 
appears that, in N = Z nuclei, one obtains sufficiently good results by keeping only the 
isoscaler fields. For T / states, the isospin projection is definitely needed, and results 
obtained with this procedure will be presented elsewhere. 

As an example of realistic shell model calculations, we first consider 24 Mg with the USD 
interaction ]nj. Figure [I] shows energies and expectation values of J ■ J, for six low-lying 
states as a function of the QMCD basis dimension compared with the exact values. In this 
case we start with five significant fields, and eventually all 144 T = one-body operators 
are activated. In the process of J-drive, three values for 9 y are employed. For 800 QMCD 
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bases, the ground state energy becomes —86.91 MeV, while the exact value is —87.08 MeV. 
The dimension of the m-scheme shell-model basis for the ground state is 28 503. Thus the 
number of bases is reduced by a factor 1/35 with a loss of accuracy of only 0.17 MeV in the 
total energy. The error due to the truncation of the Hilbert space (systematic error) does 
not exceed two hundred keV in the ground-state energy in the present calculations. 

Table | shows the lowest three energy levels, where the QMCD results for 100, 400 and 
800 basis dimensions are listed together with the exact results. One finds a remarkable 
agreement between the QMCD and exact values. Note that the accuracy of these excitation 
energies is better than that of absolute energies. In fact the deviations are less than 0.15MeV 
with only 400 basis states. In the same table, several E2 transition matrix elements and 
quadrupole moments are compared with the exact values. It can be found that several 
in-band transition B(E2) values are reproduced well with only 100 basis states, and other 
matrix elements are also obtained with 400 basis states. Thus the QMCD method turns out 
to be useful, especially for the study of low-lying collective states. 

We now proceed on to full pf shell calculations. We have confirmed the feasibility of 
the QMCD method by comparing its results with the exact ones [TJ for 48 Cr with the KB3 
interaction ]TT[, as will be presented elsewhere. In this Letter, we discuss 64 Ge. The m- 



scheme dimension of the M=0 space is 1087455 228, which is the second largest one for the 
N = Z even-even pf shell nuclei. It is larger than the dimension for 48 Cr by a factor of 
about 550, and the exact diagonalization is hopeless in the near future. This nucleus is one 
of the proton-rich N = Z unstable nuclei, and experimental data suggest that it is 7-soft. 
Thus, it is quite interesting to investigate whether we can reproduce such a structure by 
using a realistic interaction, the validity of which has been examined at least for the lower 
part of the pf shell. In this Letter, we adopt the FPD6 interaction fll2"| . This interaction is 
derived by fitting experimental data in the mass range 41-49, and is suggested to be suitable 
for describing nuclei in the upper pf shell ||T3| . 



In Fig.|2|, calculated low-lying spectra are compared with experimental data. It is re- 
markable that the calculated levels show a rather good agreement with experiment without 



any adjustment. The 7-soft nature is also evident in the calculation. The calculated ratio 
of excitation energies of 2\ to 2\ is 1.9 and that of 4+ to 2f is 2.6. Experimentally these 
ratios are 1.75 and 2.27, respectively. The relative magnitudes of B(E2) values are shown 
in Fig.2. With e p = 1.33e and e n = 0.64e, B(E2;2f -> 0f)=5 x 10 2 (e 2 fm 4 ) is obtained, 
which corresponds to ~ 0.28. The B(E2) values of the 4f — > 2± and 2^ — > 2± transitions 
are about 1.3 times larger than that of B(E2;2i — > Of), suggesting 7-softness. We obtain 
B(E2;2^ -> 0t)/B(E2;2t -> 2f) ~ 2xl0~ 3 , which is quite small similarly to the experi- 
mental value, suggesting 7 ~ 30° in triaxial deformation models 0. Calculated quadrupole 
moments appear to be small (typically | Q | < lOefm 2 ), consistently with 7-softness. 

The convergence of the results in Fig.2 has been examined by several calculations with 
different stochastic parameters. The typical deviation among different calculations (statis- 
tical error) is about 100 keV for the 2^ energy level, for instance. The discrepancy between 
theoretical and experimental results comes partly from the systematic and statistical errors 
in the present method, and partly from the interaction. The former one is being reduced by 
improving the method. 

Typical occupation numbers of fj/2, P3/2, P1/2 and f$/2 orbits are 15.1, 2.6, 0.8 and 5.5, 
respectively, for low-lying states. We can see that more than six nucleons are excited from 
the (/7/2) 16 (P3/2) 8 configuration, and that even f 7 / 2 is active. One sees that all these four 
orbits are mixed. Because of the huge basis dimension mentioned before, the conventional 
shell model diagonalization is impossible. 

The QMCD method can generate, in principle, all basis states which are needed to de- 
scribe the exact eigenstate. It is free of the assumption of some specific collective coordinates 
as in the usual Generator Coordinate Method. In addition we can take into account various 
states around many different local energy minima, which is difficult in variational approaches 
with multi-Slater determinants or multi HFB states. 

In summary, it has been shown that large-scale realistic shell model calculations can be 
carried out by the QMCD method. The QMCD method has been improved considerably 
with respect to (1) the sampling of auxiliary fields based on the local energy minima, (2) 
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the selection of dominant fields, and (3) the explicit implementation of kinematic symmetry 
requirements. Several low-lying states of large systems have been described in terms of small 
numbers of QMCD basis states with the accuracy of several hundred keV in total energies. 
The accuracy of excitation energies and E2 transition matrix elements is much better. Such 
capability of describing low-lying states is the major advantage of the present method over 
the SMMC method. The present results demonstrate that the shell model calculations 
with full valence shell configurations have become feasible by the QMCD method, shading 
light upon the structure of nuclei even beyond the pf shell with more direct relation to the 
effective nucleon-nucleon interaction. The minus-sign problem seems to be absent in the 
QMCD method, and, hence, any effective two-body interaction can be used as it is. 

We acknowledge Professors B. A. Brown and A. Poves for providing relevant two-body 
matrix elements. We are grateful to Drs. W. Bentz and A. Gelberg for reading the 
manuscript. This work was supported in part by Grant-in- Aid for Scientific Research on 
Priority Areas (No. 05243102) from the Ministry of Education, Science and Culture. A 
part of this work was carried out on the VPP500 computer at RIKEN as a part of the 
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TABLES 

TABLE I. Comparison between the QMCD and the exact results for excitation energies (MeV), 
B(E2) (e 2 fm 4 ) and quadrupole moments (efm 2 ). The effective charges e p + e n = 1.78e are used 
0- 



observable 




QMCD dimension 




exact 




100 


400 


800 




E x {2+) 


1.50 


1.54 


1.53 


1.51 




4.33 


4.23 


4.18 


4.12 




4.54 


4.50 


4.46 


4.37 


B{E2;2\ -> Of) 


74.1 


73.2 


74.2 


76.1 


B(E2;2+ - 0+) 


7.1 


7.2 


7.1 


6.8 


B{E2-2+ 2+) 


12.1 


16.8 


16.2 


16.6 


B(E2;4+ 2+) 


103.8 


102.6 


102.0 


101.1 


B(E2-M -> 2+) 


1.8 


0.4 


0.5 


0.5 


Q(2i~) 


-18.7 


-18.4 


-17.9 


-17.1 


Q(2+) 


18.5 


18.4 


18.1 


17.3 


0(4) 


-21.1 


-21.5 


-21.2 


-20.8 
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FIGURES 

FIG. 1. (a) Energies and (b) expectation values of J • J of the lowest six states of 24 Mg plotted 
as a function of the QMCD basis dimension, with N t = 20 and A/3 = 0.07 (MeV -1 ). The exact 
values are shown by symbols. Different symbols indicate different angular momenta. 

FIG. 2. Experimental and calculated energy levels of 64 Ge. The QMCD parameters are Nt = 40 
and A/3 = 0.06 (MeV -1 ). The arrows designate E2 transitions with 5(£'2)'s indicated by their 
widths. 
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